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Abstract 

We consider the motion of a collection of fluid loaded elastic plates, situated hor- 
izontally in an infinitely long channel. We use a new, unified approach to boundary 
value problems, introduced by A.S. Fokas in the late 1990's, and show the problem 
is equivalent to a system of 1-parameter integral equations. We give a detailed 
study of the linear problem, providing explicit solutions and well-posedness results 
in terms of standard Sobolev spaces. We show that the associated Cauchy problem 
is completely determined by a matrix, which depends solely on the mean seperation 
of the plates and the horizontal velocity of each of the driving fluids. This matrix 
corresponds to the infinitesimal generator of the Co-semigroup for the evolution 
equations in Fourier space. By analysing the properties of this matrix, we classify 
necessary and sufficient conditions for which the problem is asymptotically stable. 



* a.c.l.ashton@damtp.cam.ac.uk 



1 Introduction 



In this paper we are concerned with the motion of a collection of elastic plates residing 
in a horizontal channel. Between each plate is a fluid with mean flow in the horizonal 
direction which drive the motion of the elastic plates. 

The physical significance of the problem is discussed in detail in [7] , where the authors 
give both a theoretical and experimental study of the flapping modes of two elastic bodies 
of finite length. In this paper, following jj, we do not impose any external length scale. 
As such, the results are particularly applicable to cases in which the elastic bodies are 
long: an example being the motion of underwater cables. 

At the linear level, many approaches to similar problems see authors choosing to work 
with the case in which the channel is infinite in both the vertical and horizontal direc- 
tions. This approach often serves as a good approximation to some physically interesting 
scenarios, and also simplifies the analysis if working with a combination of Laplace and 
Fourier transforms. The use of the Laplace transform is generally inappropriate for the 
following reason: (a) the plates are assumed thin and as such disturbances can travel at 
arbitrarily high speeds and (b) one cannot immediately rule out the existence of instabil- 
ities that grow exponentially in time. In this paper, a new unified approach to boundary 
value problems is employed [6]. This approach is direct, and removes the need to invoke 
the Laplace transform. 

We choose to work with the more physically relevant case in which the vertical extent 
of the channel is finite and avoid the use of the Laplace tranform. The problem is 
reformulated in terms of a spectral parameter k, and in the spectral space the evolution 
equations for the amplitutdes of the plates take a simplified form. 

The problem we study is non-dimensionalised, with the only remaining parameter 
being the speed of the horizontal flow between each plate. Referring to [I] for details of 
the re-scaling scheme, the dimensionaless parameter U arises in the form: 

Here -B* is the bending stiffness of the plate, to* is the mass per unit length of the plate, 
is the density of the fluid and {/* is the horizontal velocity of the flow. In practice, the 
parameter U will be small and has served as a basis for asymptotic approaches to the 
problem in many papers. For instance, in 8 the author gives the following example: if 
the plate is made of 2cm steel, and the fluid is water, moving at around IOtos -1 (which 
can serve as an approximate upper limit), then U ~ 0.05 <C 1. 

We first work on a formal level to derive explicit solutions to the underlying problem, 
then prove what properties these solutions have. Frequently we will refer to the standard 
Sobolev space H^(R m ) = W^ 2 (R m ), where: 

^(R m ) = {d a f € L p du (R m ),0 < \a\ < k}. 

and L p dv (Yl m ) is the usual space of Lebesgue integrable functions on R m with norm 

i/p 



I/I p d^ 

The space of k times continuously diffcrcntiable functions from R m to R will be denoted 
C (R m ), and the Schwartz space iS(R m ) is defined to consist of smooth functions / from 
R m to R such that sup^. \x a d^f \ < oo for all multi-indices a, (3. 



2 




j r i+1 

-JV(r«+i) 



ri 

U 

■ ■ — B~ 

Figure 1: The geometry of the problem. 

1.1 The Governing Equations 

The problem we study concerns the motion of a collection of n clastic plates, lying hor- 
izontally, each driven by a mean flow with horizontal velocity Ui . We assume each plate 
can be described by a single valued function r]i(x,i) and we denote the corresponding 
surface by r, ; : 

r< = {(x,y) G R 2 : y = m(x,t)}, t > 0. (1.1) 

with upward normal N(Ti) = (—d x r)i, 1). The plates are situated in a horiztonal channel 
O, bound above and below by the surfacaes B + and B~ : 

B ± = {(x,y) GR 2 : y = ±h a + h ± (x)}, (1.2) 

where h > and h± G 5(R). We define Oj, 1 < i < n — 1, to be the open, connected 
region bounds between the surfaces I\ and so that: 

= {( x ,y) G R 2 : f]i{x,t) <y < r)i + i(x,t)}, t>0. (1.3) 

The connectedness condition means that for i > j, we have r}i(x,t) > rjj(x,t). In 
addition, we define the distinguished regions: 

n = {(x,y) eR 2 : -h + h-(x) <y <<q x {x,t)}, t > 0, (1.4) 
ft„ = {(x,y) GR 2 :r) n (x,t) < y < +h + h+(x)}, t>0. (1.5) 

With these definitions it is clear that = U™ =o ri,;. A diagram for the geometry of the 
problem is given in Figure |1.1| 

Each fij is occupied by an incompressible, irrotational fluid travelling with mean 
velocity Ui in x-direction. We let 4>i{x,y,t) denote the velocity potential for the pertur- 
bation from the mean flow in f2j, and as such: 

Acf)i = in f2i, (1.6) 
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which follows from the incompressibility condition. The evolution of the surface I\ is 
determined by beam equation: 

dtVi + dtVi =Pi-i~ Pi on Ti, (1.7) 

where pi = pi(x,y,t) is the pressure of the fluid in fij. In addition, on the upper and 
lower sides of Ti we have the Bernoulli conditions: 

» ^tMl^l = ~~* ) onT, (1.8) 
where we have introduced the differential operator djj = dt + Ud x for convenience of 



notation. Using the pair of Bernoulli conditions in (1.8), we can eliminate the pressure 
terms in ( |1.7| to find the single boundary condition on the surface I\ in terms of the 
functions {r/i, fa-i}: 

d1 m + 8% + du^fo-i - du^i + §||Wi_i|| 2 - §||V<&|| 2 = on IV (1.9) 

We also have two kinematic conditions on each Ti, since fluid particles initially on 
the surface must remain on the surface. 



(Ui-i + d x (f>i-i,d y <j)i-i) ■ N(Ti) = d t ra 
(iUi-i + d^i^dyfa^-iUi + d^dy^yNiTi) = 



on I\ (1.10) 



To summarise, the nonlinear hydrodynamic problem studied in this paper is governed 
by the following n + 1 dynamic boundary value problems: 

The internal domains Qj, 1 < i < n — 1: 



A0 4 = 


= 


in tti, 


(1.11a) 


(U i + d x cf> i ,d y ct> i )-N(T i ) = 


= d t r)i 


on Ti, 


(1.11b) 


d 2 tVl + dtm + du^i-r + i||W>i-i|| 2 = 




on Ti, 


(1.11c) 


{U i + d x <t>ud v <f>i)-N{T i+1 ) = 


d t r) i+1 


on r i+ i, 


(l.lld) 


d? Vl+1 + d i xVl+1 + d (A cb. l + lliv^ii 2 ^ 


: % l+1 ^+i + 5l|V^ 1+ i|| 2 


on T i+1 . 


(l.lle) 


The fact that two conditions are given on 


the upper and lower parts of d£li is 


a conse- 



quence of the fact that the surfaces Ti are unknown. 
The top domain fi„: 





A0„ 


= 


in fl n , 


(1.12a) 




-d x <t> n ,d y <t> n )-N(r n ) 


= <9 t 77„ 


on r„, 


(1.12b) 


dtVn + dtvn + d Un 


_!^n-l + i||V0„_i|| 2 


= d Un <f>i 


+ i||V0„|| 2 on r„, 


(1.12c) 


(Ui- 


f d x <t>i,dy<t>i) -N{B + ) 


= 


on B+. 


(1.12d) 



Here only one boundary condition is given on B~ , since this boundary is fixed and 
determined in terms of the given funciton h + (x). This condition is the usual Neumann 
boundary condition reflecting the fact that the fluid cannot escape through B~ 
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The bottom domain fl<j: 







= 


in fi , 


(1.13a) 


(U + d x (j>Q,dy 


M • N(Tx) 


= dtrj! 


on r l5 


(1.13b) 


9tVi + d*r}i + d Uo <j> 


+ fllWof 


= d Vl 0i H 


-fHWif ottTi, 


(1.13c) 


(U + d x (/>o,d y 4 


'o)-iV(B-) 


= 


on B~. 


(1.13d) 



Again, only one boundary condition (Neumann) is given on the fixed surface B , since 
this surface is given in terms of the known funciton /i_ (x) . 

In addition to each system of equations, we require that the amplitudes decay to 
zero at large x, and similarly we assume cj> and its derivatives decay to zero for large x, 
uniformly in y. With this assumption the "effective" part of each <9f^ is I\ (bottom) 
and IV|_i (top), since all the functions are zero on dili \ (Ti U r^+i). 



1.2 The New Dependent Coordinates 



In each of the boundary value problems listed above, it is clear that whatever the har- 
monic function (fn is, it will be determined by the values of it and its derivatives on the 
boundary <9fV It is convenient then, to reformulate the problem in terms of the the 
potentials evaluated on IV We do this by introducing the new functions ^ (x, t): 



£f(x,t) = <j>i\r i+1 , 
£~(x,t) = <f>i\ Fi , 



(1.14a) 
(1.14b) 



where </>|r; = <fr(x,r]i(x,t),t) etc. So £ + (x,i) is the potential evaluated on the top of 
the domain f2j (i.e. I\ + i), and £~(x,i) is the potential evaluated on the bottom of the 
domain (i.e. Ti). The chain rule gives rise to the following on I\ + i: 



9xt£ = (d x (l>i 



- d x T] i+1 d v (pi)\r i+1 , 
■ d t r] l+1 d y (f>i) |r i+1 , 



and similarly on IV 



d x S,i = (9x4>i + d x Tjidy4>i)\vi 
d t £~ = (dtfa + d t f]id y (j)i)\r z . 



(1.15a) 
(1.15b) 



(1.16a) 
(1.16b) 



It is possible to express all the derivatives of 4>i on the boundary d£li in terms of these 
new functions. For example, (1.10) and (1.16 1 yield the following non-singular set of 
equations for the derivatives of < 



on r, : 



d x r]i 


1 


0" 




~dx$i 




'duiVi 





dtVi 


1 




d y 4>i 






1 


dxVi 


0_ 




dt<t>i- 




AC. 



A similar set of equations hold for the derivatives of 1 
and using the notation (X) = (1 + X 2 ) 1 / 2 gives: 



on IV (1.17) 
on Ti + i. Solving these equations 



~dx4>i 


1 


dy(j)i 


(dxVi) 2 



dx&dxVi + du t Vi 
{d x 7]i) 2 dti~ - dtrjiidu^i + d x T]id x £~). 



on Ti 



(1.18) 



5 



and 



d x (j>i 
d y (j>i 



(d x r). 



■i+U 



dxif - dxVi+iduiVi+i 



on r lH 



{d x r]i + i) 2 d t £,t - dtm+xiduim+i + d x r) i+1 d x £, + ) 

(1.19) 

These expressions indicate that it is enough to know the functions £f(x,t) to reproduce 
the potential 4>i(x, y, t) on f2j. The problem now consists of determining each amplitutdc 
rji(x,t) and the potentials each dSli. The total number of unknowns is: 



(# of amplitudes) + (# of potentials on each I\) + (potentials on S d 



3n 



Our aim is now to reduce our complex boundary value problem into a system of 3n + 
2 equations for the unknowns. The Bernoulli condition (1.9 1 on each I\ will give n 
equations in terms of the and {rji\, so we are left to produce 2{n+ 1) independent 

equations, i.e. two equations from each f2 i; so that the problem is well posed. 



2 The Non-Local Formulation 

In this section we present a non-local formulation of the problem that involves translating 
the harmonicity of each (f>i, along with the dynamic (or Neumann) boundary conditions 
on d£li, into a pair of integro-differential equations in ^ and r\i. This is achieved by 
constructing a so-called global relation El H] ■ The global relation is a consequence of 
the following simple lemma. 



Lemma 1. Suppose that the functions u and v are harmonic in Q C R 2 . 
following holds: 



d x (d y u d x v + d y v d x u) + d y {d y u d y v — d x u d x v) = 



Then the 



(2.1) 



for each (x, y) E SI. 



This can be verified by expanding out the left hand side of (2.1 1 to find: 

(dyv)Au + (d y u)Av 
which vanishes in 17 if both u and v are harmonic in f2. 



In what follows we shall be integrating expressions such as (2.1l over each fli, and 
as such we need conditions on the functions {r]i,4>i} and their derivatives to decay suf- 
ficiently rapidly as |x| — > oo. At this stage we make the assumption, a priori, that the 
functions have the necessary conditions so that the integrals converge. It is possible to 
make these conditions precise using standard Sobolev estimates, but it is preferable at 
this stage to work on a formal level. 

Suppose fceR and define v g C°°(R 2 ) by: 



v(x, y) = exp(— ifcx + ny), 



±k. 



Clearly v is harmonic and bounded in each f2j . It follows from an application of lemma 
[TJthat the following holds in each fi,: 

dx (e-ikx+ny^ dx(f)i _ ikd^.fj + d y (e- ikx+ ^{Kdy^+ikd x ^)) = 0- (2.2) 

We now look at internal, top and bottom domains in SI seperately. 



G 



The internal domain fij, 1 < i < n — 1: 



We integrate (2.2 1 about an internal domain f2j. An application of the divergence theo- 
rem gives: 



e- ikx+Ky {Kd x <f>i - ikdyfr, Kdyfa + ikd x <t>i) ■ N(Ti) dx 



-ikx+w^fc _ ifc^^, ndy^ + ikdxfa) ■ N(T l+1 ) dx = 0, (2.3) 



where we have discarded the contributions from \x\ — > oo, since the fields arc assumed 



to vanish there. The contribution in (2.3 1 from I\ is given by: 



-ikx+ KVt ^ K (q^. _ d x Tfid x (j>i) \ Ti + ik {dxfa + d x r]id y (j) t ) | r< 



da;. 



We see that the first term in the integrand is du^i by (1.11b), and the second is 9^ 

(2.4) 



by (1.16). The expression becomes: 

Performing a similar calculation, we find the contribution from to be: 

- J e- ikx+Kr >*+i (ndu^i+i + ifca^+) dx. (2.5) 



Combining equations (2.4 1, (2.5 1 and (2.3 1, we find the following 1-parameter family of 



integrodifferential equations: 

-ite^ww ( KduiTj . + ikd^l) - e^ 1 (ndu.rii+t + ikd x &) ) dx = 0. (2.6) 



Recalling that k — ±fc, (2.6) gives two equations, which can be added and subtracted to 
give the following result. 



Proposition 1. The incomplete boundary value problem in (l.llal, (1.1 lb I and (1.1 Id I 

is equivalent to the the pair of integro- differential equations: 



~ lhx (dutVi sinh(fcr?j) + id x ^ cosh(fc?7 4 ) 

- duiVi+i smh(kr] i+ i) - id x £+ cosh(kr) l+1 )jdx = 0, (2.7a) 



and 



J e lkx (dujh coshikrit) + id x ^ sinh(fc?7j) 

- d Ut ri i+ i cosh(fcr] i+ i) - id x £f sinh(fcry 4+ i)^da; = 0, (2.7b) 

which are valid for k 6 R. 

Remark 1. We note that if both 7^ + i and r\i are small, in an appropriate sense, then 



(2.7a I and the completeness of the Fourier transform yield d x ^ + ~ 3 X £ . This is consis- 



tent with what one would expect, because over short length scales the speed of the flow 
in the x-direction would remain approximately constant in the absense of viscosity. 
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The top domain f2 ra : 



Integrating (2.2 1 over f2„ and using the divergence theorem we find: 



In analogy with (2.4 1 we find the contribution from r„ to be: 



(2.9) 



The contribution from B + , however, takes a simpler form owing to the fact that is 



boundary is fixed. Using the Neumann condition ( 1 . 12d I with N(B + ) = (—d x h + , 1) and 



the definition £^"(a;, t) — <j>(x, ho + h + (x), t) we find the contribution from B + to be: 

-ik [ e- [kx+K( - ho+h ^d x ^dx. (2.10) 



Combining ( 2.9 ) and (2.101 we find the following 1-parameter family of integro-differential 
equations: 



-ikx ^ ( KdUnVn + i k d x £-) - ike< h ° +h +)d x &) dx = 0. 
In analogy with Proposition [l] we arrive at the following. 



Proposition 2. The incomplete boundary value problem in (1.12al, ( 1 . 12b I and (1.1 2d I 

is equivalent to the the pair of integro-differential equations: 

e -ikx fd Un r] n smh(ki] n ) + id x £~ cosh(krj n ) - id x £+ cosh(k(h + h + ))jdx = 0, (2.11a) 



and 



J e-' lkx (d Un r) n cosh(fcr?„) + id x {~ sinh(fe7 ? „) - id x £+ sinh(k(h + h + )))dx = 0, (2.11b) 
which are valid for k € R. 

Remark 2. We see that if B + is flat, i.e. h + = 0, then we can eliminate £+(cc,i) from 



equations (2.11 1. Multiplying (2.11a) by sinh(/c/i ), ( 2.11b I by cosh(kh ) and subtracting 
yields the single global relation: 



du„VnCOsh(k(r) n - ho)) + ikd x t; n sinh(fc(^„ - h )))dx = 



(2.12) 



The bottom domain f2. 



o- 



The derivation of the results for the bottom domain is analogous to the previous calcu- 
lation for the top domain. 



Proposition 3. The incomplete boundary value problem in (1.13al, ( 1 . 13b I and ( 1 . 13d I 

is equivalent to the the pair of integro-differential equations: 

j e~ ikx (d Uo r]i sinh^i) + id x $ cosh^i) - id x ^ cosh(fc(/i - h-)))dx = 0, (2.13a) 
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and 

e-' lkx (du^coshik^) + id^ sinn^) + sinh(fc(/i - h-)fjdx = 0, (2.13b) 
which are valid for k £ R. 

Remark 3. In analogy with the result in Remark[2] we note that in the case i3_ flat, we 



can eliminate £ (x,t) from equations (2.131. Multiplying (2.13a) by sinh(fcft, ); ( 2.13b I 



by cosh(fc/io) and adding yields the single global relation: 

~ ikx (du m cosh{k(rn + ho)) + ikd x $ sinh(%i + /i )))dz = (2.14) 

The non-local formulation of the problem has given two explicit equations in each 
^ij < i < n for the functions }■ Together with the n Bernoulli-type conditions, we 
have: 

2n + 2 + n = 'in + 2 

equations for the 3n + 2 unknowns, so our problem is well-posed in this sense. The 
important feature of our 3n + 2 equations is that they are explicit in the functions 
{rii,^}, which makes analysis of them much simpler. For instance, they allow direct 
use of perturbation expansions, direct application of standard Sobolev estimates and 
direct application of standard asymptotic techniques - we refer to [2] for examples of 
such techniques used in a related non-local formulation for rotational water waves. 

We end this section with the following result, which summarises the non-local for- 
mulation for the dynamic boundary value problem outlined in ( |1.11[ ), ( |1.12 1 and (1.131 



along with the n bernoulli-type condtions on the {rj}™ =1 in terms of {^f}?=o 



Proposition 4. The solution to the boundary value problem in ( 1.11| ) ; (1.12) and (1.131 



is completely determined by the 3n + 2 functions {£j }f =0 which satisfy: 

• The 2n — 2 non-local internal equations: 

~ lkx (duiVi sinh(krn) + id x ^ cosh(Zc^) 

- du.rn+i sinh(fc77i+i) - id x $~ cosh(kr) l+1 )jdx = 0, (2.15) 

-ikx fg UiT j. C osh(fcry 4 ) + id x ^ sinh^) 

- duiVi+i cosh(/cry i+ i) - id x £+ sinh(fc^ + i)^da; = 0, (2.16) 

which are valid for k £ R and 1 < i < n — 1 . 

• The 4 non-local external equations: 



e~ ik 



e~ ik 



x (d Un r) n smh(kr] n ) + id x £ n cosh(/cry„) - id x £+ cosh(fc(/i + h + ))^dx = 0, 

(2.17a) 

x (pu n nn cosh(fc?7„) + \d x (,~ sinh(fc?7 n ) - id x £+ sinh(fc(/i + K))) dx = °> 

(2.17b) 

lkx (duoVi sinh(fcr?i) + id x ^ cosh(krii) — id x £,Q cosh(fc(/i — h_))jdx = 0, 

(2.17c) 
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J e lkx (du m cosh(fc?7i) + id x $ sinh(fc?7i) + i9 x £ sinh(fc(/i - h-))jdx = 0, 

(2.17d) 

which are valid for fceR. 
• The n Bernoulli-type equations: 

df Vi + d% + du^fa-i + H|V&_i|| 2 = duifc + IllV^H 2 on Ti (2.18) 



for 1 < i < n, where an d ^t<^i|r 4 e ^ c - are given in terms of via (1.181 



and (1.191 



3 The Linear Problem and Well-Posedness 

In this section we consider the linearisation of the equations appearing Proposition [4] and 
the corresponding Cauchy problem for the linearised equations. We denote iji — hi + rji 
where the {hi} are constant and ||9 x £i Hl 00 an d H^IU 00 are small so that quadratic terms 
and beyond may be neglected. For simplicity we restrict attention to the case in which 
B are flat. It is straight forward to adjust the analysis to deal with the case in which 
the upper and lower parts of dfl vary according to a funciton in the Schwartz class. 

We use / to denote the Fourier transform of the function /, and drop the tilde from 
fji for notational convenience. It will also be useful to denote the upper and lower parts 
of the channel by y = h n +i (top) and y = ho (bottom). Discarding quadratic terms and 
higher, the linearised equations are written as: 

• The 2n — 2 internal equations: 

(dtfji + ikUifji) sinh(fc/ij) - (dtfji+i + ikUifj i+ i) sinh(fc/i i+ i) 

+ fcQ cosh(fc/ii+i) — fc£~ cosh(khi) = 0, (3.1a) 



(d t fji + ikUifii) cosh(khi) - (d t fj i+ i + ifeUj^+i) cosh(kh i+ i) 

+ k^ s,inh(kh l+1 ) - kii sinh(khi) = 0, (3.1b) 

for 1 < i < n. 

• The 2 external equations: 

(dtVn + ikU n f) n ) cosh(fc(/i n+1 - h n )) + fc|~ sinh(fe(/i n+1 - h n )) = 0, (3.2a) 
{d t f]i + ikU fio) cosh(fc(/ii - ho)) + k££ sinh(fc(/i a - ho)) = 0. (3.2b) 

• The n Bernoulli-type conditions: 

+ k% + (ddU + iAE/i-i&tx) = {di~ + ikU^), (3-3) 

for 1 < i < n. 

Note that we have eliminated and £q using two of the external equations (see Remarks 
[2] and [3]) leaving 3n equations for 3n unknowns. Equations (3.1l and (3.2) allow us 



to express the in terms of {fji}, which can then be used in (3.3 1 to determine 
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the governing equation for the {fji}. Introducing the functions tp(k,t) = (fji, . . . ,f] n ) T , 
n(k, t) — (dtfji, . . . , d t r)n) T and Ahi — h i+1 — hi, the governing equations are: 



csch(fcAfej) 
k 



7Tj+l + 

2i£/ i csch(fcA/i j ) 
kllfcsch(kAhi) 



-i . coth(kAhi) . coth(kAhi-i) 
1 n E I E 



■hi + 



csch(kAhi-i) 
k 



2iC/ 4 coth(A;A/i 4 ) + 2i£/ i _ 1 coth(fcA/i l _ 1 ) 
fc 4 - fc[/f_ 1 coth(fcA/i^i) - kUf coth(fcAZii) 



+ 



2i[/ i _iCsch(fcA/i i _i) 



fcU'i 2 _ 1 csch(fcA/ij_i) 



for 2 < i < n— 1, where the dot denotes differentiation with respect to t. These equations 
are supplemented with the extra equations n equations tp = ir. The problem becomes 
one of linear algebra: it is straightforward to see that the governing equations will take 
the form: 



Mi(k)dt 



= M 2 {k) 



(3.4) 



where Mi(k) £ Mat2„(C), i = 1,2, whose entries are meromorphic in k. As such, if the 
matrices Mi and M 2 are determined explicitly, one can immediately write down the 
solution to the initial value problem: 



= exp(tM){k) 



(3.5) 



where M. = M.^[ 1 M. 2 and the vector on the RHS contains the initial data. We see 
that the Cauchy problem for the linearised equations is completely determined by the 
matrix A4(k), which corresponds to the infinitesimal generator of the Co-semigroup for 
the evolution equation in the Fourier space. In what follows we prove several results 
regarding the matrices M\,M 2 . 



Lemma 2. The matrices Mi(k) and M 2 (k) take the form 
^iW = (o"^)) and M 2 (k) 



i„ 

fc£(fc)-fc 4 I„ -2iC(fc) 



where A(k),B(k),C(k) £ Mat„(R) are symmetric, tridiagonal, meromorphic in k and 
real for k e R. The matrix A(k) is given by: 







fc+coth(fcA/ti)+coth(fcAfo ) csch(fcAZii) 
k k 
csch(fcAfei) fc+coth(fcAfe 2 )+coth(fcA/ti) csch(fcA/t 2 ) 

k k k 

csch(fcAfe 2 ) fc+coth(fcAfo 3 )+coth(fcAfe 2 ) 
~ k k 







The matrix B(k) is given by: 

Ug coth(kAh )+U^ coth(feA/ii) —Uf csch(fcA/i!) 

-(7 1 2 csch(feA/ii) Ul coth(feA/ii)+(/2 coth(feA/i 2 ) -t/f csch(fcAh 2 ) 

-!7 2 csch(fcA/i 2 ) Ul coth(fcA/i 2 ) + (7| coth(fcAfi 3 ) 



The matrix C(k) is given by: 

~ U coth(fcA/i ) + t7i coth(fcAfti) -!7icsch(fcA7ii) 

-!7icsch(feA/ii) Ui coth(kAh 1 ) + U 2 coth(fcA?i 2 ) -C/ 2 csch(fcA/i 2 ) 

-U 2 csch(kAh 2 ) U 2 coth(kAh 2 )+U 3 coth(feAft, 3 ) 
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Whilst the general form of M.(k) is complicated, many of the entries become expo- 
nentially small beyond \kh\ > 0(1), where h = max^ Ahi. This makes the asymptotic 
analysis of M(k) extremely simple. In dropping the exponentially small terms, the Ahi 
are removed from the problem and we deduce that the short wavelength dynamics are in- 
dependent of the distribution of the elastic plates within f2. This is expected on physical 
grounds: the separation will effect modes with wavcslcngths of comparable size. 

Lemma 3. The determinant of M.\(k) is non-zero for each k G R. 

Proof. It is enough to prove that A(k) is invertible, and we do this by proving „4 ^ 0, i.e. 
A is positive definite. Denoting the diagonal entries {a,i}™ =1 , and the sub/sup-diagonals 
by {-bi}?~i, we have for X £ R": 

n n—1 

X ■ (AX) = aiXf - 2 biXiX i+1 

i=l i=l 

n-1 

> (at - b{)Xl + (a n - b n ^)X 2 n + ^(a, - b, - b^Xf, 

which follows from AM-GM inquality. It is straight forward to prove that each of the 
coefficients of the Xf are strictly positive, i.e. 

\(k + coth(fcA/ij) - csch(fcA/ij) + coth(fcA/i i+ i) - csch(fcA/i l+ i)) > 

so we conclude ^1^0. From this we deduce that det(A^i) > for k G R. ■ 

An immediate corollary to Lemma[3]is that M\(k) is invertible, and what is more, 
A4i 1 (k) is analytic on the real line. This follows from the previous lemma and the fact 
that the entries of A(k) are analytic on the real line. 

Lemma 4. The matrix M(k) is a analytic function of k in a neighbourhood of the real 
line whose spectrum becomes purely imaginary as \k\ — > oo. 

Proof. Our previous discussion establishes the analyticity of Aii 1 (k) on the real line. It 
is clear that M 2 (k) is analytic away from the origin, and a straightforward computation 
of the relevant Laurent series about k — establishes analyticity of (M. 1 ~ 1 Ai2)(k) at 
k = 0, so the first claim is established. For the second, some standard asymptotic 
expansions for large, positive k reveals: 

A(k) ~ I„, (3.6a) 
B(k) ~ diag(t/ 2 + Ul . . . , + Ul), (3.6b) 
C(k) ~ diag(C/ + Ut,..., U n -i + U n ). (3.6c) 

This gives us the asymptotic form of the characteristic polynomial: 

n-l 

det(M(k) - AI 2 „) ~ (-1)" l[ (A 2 + 2i(U m + U m+X ) - k(U 2 m + U 2 m+1 ) + fc 4 ) . 

m— 

Since the characteristic polynomial and hence the eigenvalues of M.(k) depend continu- 
ously on the elements of M.(k), we deduce that the eigenvalues become purely imaginary 
as k — > oo. A similar argument holds for k — > — oo. I 

We now address the issue of well-posedness: given initial data in a specified function 
space, what can be said in regards to the continuity of the map from the inital data to 
to the solution for times t > 0. 



12 



Theorem 1 (Local well-posedness in L^ a ,(R) x2n ). Let the n elastic plates be distributed 
throughout the channel fl. If the initial amplitudes rjo and velocities dtr]o belong to 
H^ X (R), where s > 2, the Cauchy problem on R x (0,T) is locally well-posed in L\ x {Yl), 
so that: 

\\(r),d t r))\\L*<C n (T)\\(r, ,d tn o)\\ H * 
where C n (T) is a constant depending on the time domain of the problem and the number 
of elastic plates. We have used the norm: 



n 



2 

i\\L- 



and similarly for \\ ■ \\h s ■ 

Proof. By Parseval's theorem it is enough to prove the relevant inequality in the Fourier 
space. We proceed to make a bound on exp(A4t)(k). We have already proved in Propo- 
sition [4] that M.{k) is analytic on the real line, so we only consider the large k behaviour 



of exp(A4t)(k). From equation (3.6 1 we see that M. can be written as: 

T(k)+£(k) 

where V(k) has purely imaginary spectrum and \\£ (k)\\ is arbitrarily small for k suffi- 
ciently larg^] Now recall the following estimate from semigroup theory: if || exp(7 : 't)|| < 
ce^ and M = V + £ , then: 

\\exp{Mt)\\ < ce^ +c " £ ll)*. (3.7) 
for some c > 0. A straightforward computation reveals: 

|| exp(Vt)\\ < nk 2 Vt. 



Choosing |fc| large enough so that ||£|| < e, the bound in (3.7 1 gives: 

||exp(Mt)(fc)|| < cnk 2 ^/te cet 

< cn(k) 2 ^e ctT . 

for k sufficiently large. Finally then: 

||?7|| L 2 = ||exp(tA^)(fc)ry |U= 
^cnVTe^Uk) 2 ^^ 
<C n (T)\\ Vo \\ Hs . 

where s > 2. This is the desired result. ■ 

This result is, of course, local. One could not hope for a global well-posedness for 
arbitrary mean flow velocities U, based purely on physical grounds. If, for instance, 
the magnitudes of the velocities were too great, or were in opposing directions, small 
instabilities would rapidly develop. See Figures 2 and 3 for an example of this phenomena. 

Using standard Sobolev theory, we observe the continuous embedding: 

where Cg (R) denotes the space of continuously differentiable functions which decay at 
infinity. As such, we deduce that local well-posedness is achieved if the initial data is 
continuously differentiable with decay at infinity, possibly after redefinition on a set of 
measure zero. 



1 Thc norm || • || can be any induced norm on Mat2n(C). 
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Figure 2: Numerical implementation of of solution to the linear Cauchy problem with n — 6. 
Four elastic plates are given Gaussian initial data, the seperation is uniform and the mean flow 
between each plate is U = 0.1. The plots are evaluated at t = 0,1,2,3 respectively. Note the 
amplitude of the disturbances decay, indicating short time stability of this flow configuration. 





Figure 3: The same numerical implementation as in Figure [3] but with the mean flows given 
by U = (+0.4, +0.4, +0.4, 0, —0.4, —0.4, —0.4). In contrast to the previous case, we see the 
amplitutde of the disturbances increasing indicating the instability of the flow configuration. 
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4 Stability 



The long time stability is determined by the Green's function of the underlying system 
of PDEs. Our non-local formulation has already provided us with the Green's function: 



from equation (3.5 1 we see that it is given by G £ Mat2«(C), where: 



G(x,t) = i- f exp{tg){k)dk (4.1) 

where we have defined the matrix Q £ Mat2 n (C) by: 

g(x,t;k)=ik(~)l 2n + M{k). (4.2) 

The problem is asymptotically stable if and only if G(x, t) — ► as t — > oo for each 
fixed V = x/t. It is clear that the asymptotic behaviour of G(x,t) is determined by the 
spectrum of Q, denoted a(G). In particular, the leading order term in the asymptotic 
expansion will be determined by the eigenvalue with maximal real part. If 5RA > 
for any A £ cr(Q), then the leading order behaviour of G(x, t) has exponential growth 
and the problem is asymptotically unstable. In the special case where o~(Q) £ iR, the 
leading order of G can be computed via the method of stationary phase, which gives 
G(x,t) ~ i" 1 / 2 as t — > oo, indicating a stable configuration^] Alternatively one can view 
this as a consequence of the Riemann-Lebesgue lemma. 

Theorem 2. The linear problem is asymptotically stable if and only if there is no mean 
flow, i. e. Ui = for i = 0, . . . , n. 

Proof. The problem is certainly unstable if there is at least one clement of o-{Q) with 
positive real part, so it is sufficient to prove the existence of such an eigenvalue to prove 
the problem is unstable. We see from lemma|2]that M{k) has the form: 

In 

k {A- X B) (k) - k+A-^k) -2i{A- l C){k) 



Since (A~ 1 C)(k) £ Mat„(R), we see that 9t(TrM) = and so from (4.2 1 



3t(TrG) = J2 ^ A = - 

It follows that either (a) all eigenvalues of Q are purely imaginary for each k £ R or (b) 
the real parts of at least two members of cr(Q) have opposite sign for at least one value 
of k £ R. If (b) is true, we are done, since there must be at least one member of <r(Q) 
which has positive real part. 

Now consider the case in which (a) holds, so that a(Q) C iR, or equivalently a(A4) C 
iR. For A £ a(A4), set A = i/x, where /j, £ R. A straightforward computation gives the 
characteristic polynomial for M.: 

XM (i M ) = det (n 2 A + 2fiC + kB- fc 4 I„) . (4.3) 



If we show that the matrix in parenthesis appearing in (4.3 1 is positive definite for some 



range of k (for all real fi), then |xai(im)I > in that range, for all real /i, and so there 



2 Here we have assumed the relevant matrix is diagonalisable, which turns out to be the case when 
<j(6) C iR. 
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are no purely imaginary eigenvalues. Using a similar argument to that in the proof of 
lemma [3] one finds the matrix appearing in ( |4.3[ ) is positive definite for all /i e R if: 

[t/jtanh (ifcA/ij) + {7i_itanh (|fcA^_i) ] 2 

< [-/c 3 + C/ftanh(ifcA^) + C/f_ 1 tanh(|fcA/i i _ 1 )] 

x [fc + tanh (ifcA/ij) + tanh (§fcA/ij_i) ] (4.4) 

for each i = 1, ...,n. This comes from examining the discriminant of the relevant 



quadratic in /i. There is always some non-empty fc-interval in which (4.4 1 is satisfied 
(choose k sufficiently small), unless Ui = for i = 0, . . . , n, in which case (4.4) cannot 
be satisfied for any k. We deduce that a(Ad) always contains an element with positive 
real part, unless Ui = for i = 0, . . . , n in which case: 

which is well defined since A >~ 0, by lemma [3] In addition, we note that A is a Jacobi 
matrix so has n distinct eigenvalues. It follows that in the case [7^ = for i = 0, . . . , n, 
Q(k) has 2n purely imaginary, distinct eigenvalues. In particular, Q{k) is diagonalizablc 
for each k e R. We conclude that the problem is asymptotically stable if and only if 
there is no flow, in which case G(x,t) ~ i -1 / 2 as t — > oo, with V = x/i fixed. ■ 

This result can be reconciled with the well-posedness result in theorem [l] In the case 
that the U are not all zero, we know that o~(Ai) is not purely imaginary for all k £ R, 
but rather a(Ad) is becomes purely imaginary for \k\ large. As such, the optimal bound 



on || exjp(Adt)\\ was of the form in (3.7l, which led to the local well-posedness result. If 
we consider the case Ui = for i = 0, . . . , n, then a(M) C iR, and the optimal bound 
in this case is found by using: 

||exp(Mi)|| < n{M)e a( - M) \ (4.5) 
where a(M) is the spectral abscissa of the matrix Ai, defined by: 

a(M) = sup 3?A 

and k(A4) = ||V||||V _1 || is the condition number of A4, where V is the matrix of eigen- 
vectors of M. After some standard asymptotic estimates, one can show: 

ifc 2 I?i ifc 2 Iji 

V{k) ~ 

Using this, it is straightforward to show that k(M) is 0(k 2 ) for large k. Finally, using 



the estimate in (4.5 1 and the fact that a (Ad) — 0, we have: 

\\eMMt)(k)\\ < C n (k) 2 

where C n is a constant depending only on n. In summary, for that case of no mean flow, 
so that Ui = for i = 0, . . . , n, we have following global well-posedness result. 

Theorem 3 (Global well-posedness in L 2 lx (TV) x2n in the case of zero mean flow). Let 
the n elastic plates be distributed throughout the channel £1, and suppose there is no mean 
flow so that Ui = for i = 0, . . . , n. Then if the initial amplitudes r/o and velocities dtrjo 
belong to H^ x (Il), with s > 2, the Cauchy problem is globally well-posed in Lj a (R). 



1G 



Let us now examine the result in Theorem [2] We see that the condition that %ta(Q) = 
is only violated for small k, i.e. only modes corresponding to long wavelengths are 



unstable. Let K denote the interval for which (4.4 1 is satisfied, so that for k e K we 
know dia(Ai) ^ 0. To understand the instabilities it is enough to know A4 for k 6 K, 
i.e. k small. It is straightforward, although algebraically intensive, to show that the 
leading order terms in the real part of the eigenvalues of M. are of the form: 

£(U,Ah) 



mi 



-k + 0(k 2 ). (4.6) 



where C is a bilinear function of U = (Uq, U\, . . .) and Ah = (Aho, Ahi, . . .) and |fi| = 
h n+ i — Iiq is the channel width. For large seperations Ah, and mean flows with small 
velocities, this quantity is exceedingly small. As such, the growth rate of the amplitudes 
of the corresponding modes will be very slow. One would expect to have to wait for long 
times before such instabilities were realised, as indicated in Figure [3] 

We have shown that the linear problem is asymptotically stable if and only if there is 
no mean flow, however it is unlikely that in practice this stability will be realised. It is 
well known that the naive treatment of stability by means of spectra can be misleading, 
the more relevant quantity being the pseudospectrum: 

a t (M) = {A G C : \\(M - AI) — 1 1| > e" 1 }. 

The pseudospectrum contains information about how the spectrum of the associated op- 
erator changes under small perturbations, in the sense that a e (M) contains information 
about the spectrum of the operator M. + £, where \\£ || < e. See figure [I] for some plots 
of a e (M) for given e. Recall that the pseudospectrum of a normal operator Af is: 

U B *w, 

i.e. one draws balls of radius e about each element of the spectrum of AT. In this sense, 
the deviation from this form of pseudospectrum is a measure of non-normality. Sec 9 
for several different, but equivalent definitions of pseudospectra and non-normality. It 
is clear from Figure that for k — 1 the operator is approximately normal, and as |fc| 
decreases the non-normality of M. increases. Indeed, the width a e (Ai) is approximately 
8e in the final output. This indicates that whilst theoretically the instabilities of the 
long-wavelength modes grow very slowly, in practice the instabilities will grow at a much 
faster rate since 3ta e (Ai) approaches 0(1) for small k. 

At this stage we recall a result of Chrighton & Oswell [3], who proved that in the 
case of a single elastic plate driven by one mean flow (above) the configuration was 
convectively stable for U < U c = 0.074. However, the configuration in their work differs 
considerably from that considered here - most importantly the flow is not bound above 
(i.e. their channel is infinite in width). In our problem, the fluid in £li transfers vertical 
momentum betwen r\i and and in particular the rigid boundaries of the channel 

reflect back waves carrying vertical momentum. In the case of non-zero mean flow, as the 
fluid interacts with the plates it is able increase the vertical component of the momentum 
- this process continues and instabilities are liable to grow in time. However, in [3] the 
infinite vertical extend of the channel means that disturbances can "escape" from the 
channel, since the fluid can carry vertical momentum away from the plate and out to 
infinity, and the aforementioned incease in vertical momentum does not occur. One can 



hcuristically reconcile their work with ours by taking |f2| — * oo in (4.6) 
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Figure 4: Plots of the pseudospectrum of M(k) with n = 4, for k decreasing from 1 to 0. 
Plate seperations are all 10, and U = (+0.15, +0.3, 0, —0.3, —0.15). The plots are given by 
\\(M - AI)" 1 !! = e -1 , where e = 10" \ 10~ 15 , 10~ 2 , 10~ 2 ' 5 . The width of a,(M) increases as k 
decreases, indicating the non-normality of the matrix M for small k. 



5 Variations of the Problem 

In this section we briefly discuss some variations and extensions of our problem. For the 
sake of brevity, some of the details are omitted, but all the would-be calculations are 
similar to those seen in the prcqual. 

5.1 An Infinite Channel 

In some certain circumstances it is more appropriate to treat the surfaces B± at infinity. 
In this case, the analysis changes on the top and bottom domains. Recall that in the 



non-local formulation on the top domain, the initial step (2.8) involved the integral: 

e- ikx+Ky {nd x <f>i - ikdy^Kdyfa + ikd^i) ■ N(T n ) dx 

r„ 

- / e-' lkx+K y{Kd x & - ikdyfa, Kd y & + ikd^i) ■ N(B+) dx = 0. 

JB+ 

In the case of the infinite channel, we are required to take the limit y — > oo in the second 
integral and as such we need to impose k — k < 0. In this case, the integral equation in 
the top domain fi„ becomes: 



k X+ k Vn + dx = Q; (5 1} 

with the restriction k < 0. Alternatively, changing k t— > —k, (5.1 1 becomes: 

(Q UnVn + idx ^ dx = 0j k > 0. (5.2) 
The corresponding equation for the bottom domain f2 i s: 

,-ikx+km (g Uom + idx ^ dx = ^ k > o. (5.3) 

The analysis on the internal domains will be exactly the same as in the case of the finite 



channel. The restrictions that k > in equations (5.2 1 and (5.3 1 does not cause any 
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restrictions, at least in the linear case, because a real function can be recovered from the 
knowledge of its Fourier transform on k > 0. This follows from the observation: 



e lkx f(k) dk 



Akx 



f(k) 



c f(k) 



-ikx 



/(-*) 



dk 



^—ik 



x f(k) 



dk, 



which is a consequence of the reality of / and the definition of the Fourier transform. 
Since the functions in our problem are manifestly real, we see that the restriction that 
k > on the top and bottom domains is permissable. A similar restriction was made in 
[3] in which the authors prove well-posedness results for the initial value problem for a 
single fluid loaded clastic plate on both the full and half lines. 



5.2 Non-zero Vorticity 

In the presence of non-zero vorticity, one cannot invoke a scalar velocity potential, so 
the methods outlined thus far are inappropriate. In this case, one refers to the stream 
function ip defined so that ip v represents the horizontal velocity, and — ip x represents the 
vertical velocity of the fluid. For a given vorticity uj, the stream function satisfies: 

-Ai/j = lo. (5.4) 

In [2] the authors derive a non-local formulation for rotational water waves with constant 
vorticity lj — 7, by introducing the pseudo-potential tp defined as the harmonic conjugate 
to the function: 



which is harmonic by virtue of equation (5.4 1. The problem can be stated in terms of 



the harmonic function tp, and a non-local formulation follows in a similar manner as that 
seen in SI. 



5.3 Higher Dimensions 

An important generalisation of the problem is to extend it to higher dimensions, three 
spatial dimensions being the most physically relevant case. We analyse the problem in 
m + 1 spatial dimensions, with m > 1, using the coordinates: 

(xi, x m , y) G R m+1 

where y corresponds to the vertical direction. The elastic plates are now described by 
the m dimensional hypersurfaces: 

T l = {(x,y)eR m+1 :y = Vt (x,t)}, 

where x = (x\, . . . , x m ) is the "horizontal" coordinate. The mean flow in each fij is 
horizontal and can be represented by (Ui,0) G R m+1 . The scalar differential operator 
djj needs a suitable modifcation in the higher dimensional setting - it becomes: 

d v = d t + U ■ V x 

where = (d Xl , . . . , d Xm , 0) is the horizonal gradient operator. In addition, the beam 
equation for each elastic plate gets modified so that: 

dfi]i + A 2 x r). L = pi_! - pi on I 1 ,, 
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where is the laplacian in the horizontal coordinates. The relevant kinematic boundary 
conditions are modified in an obvious fashion and by introducing analogues of the 
one can express the derivatives of 4» in terms of 77 and £. Indeed, the {\7 x 4>i} satisfy the 
following non-singular equation: 

VcC," - ((>i //, iV, ;/, = (I + V x r)i <g> V x r)i)V x (()i on I\, (5.5) 
which comes from eliminating d y <pi from the higher dimensional analogues of equations 



(1.16a I and (1.10 1. From this we can express dt<f>i and d y (j)i in terms of ^ and rji on the 
surface IV The higher dimensional analogue of the identity in lemma [I] is: 

V x • ((dyu)V x v + (d y v)V x u \ + dy((d y u)(d y v) - V x u ■ V x uJ = 

which holds for all harmonic functions u, v on R m+1 . Replacing v by exp(ifc • x ± ||fc||y), 
with k S R m , the governing equations ( |1.11[ >, ( | 1 . 1 2 ) , (1.13), in the case of m + 1 spatial 
dimensions, can be reformulated in terms of a non-local problem. 



Proposition 5. The solution to the boundary value problem described by ( |1.11[ ), ( |1.12 1 
and (1.131 in the case o/m + 1 spatial dimensions is completely determined by the 3n+2 



functions {r]i}f =1 , }™ =0 which satisfy: 

• The 2n — 2 non-local internal equations: 

u '(<)> 11, sinh(||fc||7 ?l ) + ik ■ V^r cosh(||fc||^) 

- d Ui r) i+1 smh(\\k\\ri i + 1 ) - ik ■ cosh(\\k\\r) i+1 ))dx = 0, (5.6) 



Ak ' x (d^rii cosh(||fc||?7i) + ik • V x £i sinh(||fc||^) 

- 9(7 4 jy i+ icosh(||fc||j7i + i) - ik ■ siDh(\\k\\r) i+ i))dx = 0, (5.7) 



where k is the unit vector in the k-direction and dx = dxi A • • • Adx m is the volume 
form on R m . These are valid for k S R m and 1 < i < n — 1. 

• The 4 non-local external equations: 

~ ik ' x (du^nSmhiWkWrin) + ik ■ V x £" cosh(||fc||77„) 

- ik ■ V a f+cosh(||fc||(/io + h+)))dx = 0, (5.8a) 



e ik ' x (d Un r)nCOsh.(\\k\\r) n ) +ik ■ V x £ n sinh(||fc||ry„) 



ik ■ V x £+ sinh(p||(/i + h+)))dx = 0, (5.8b) 



- ifc K (% oJ?1 smh(\\k\\ m )+ik ■ V x ^+ cosh(\\k\\ m ) 

- ik ■ V x £q cosh(||fc||(/io - M) )dx = 0, (5.8c) 
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J e- ik x (d Uo rn cosh(\\k\\r n )+ik-V x ^smh(\\k\\ m ) 

+ ik ■ V x ^o sinh(||A;||(fto - h-)))dx = 0, (5.8d) 

which are valid for k G R m . 
• The n Bernoulli-type equations: 



on Ti 



(5.9) 



for 1 < i < n, where V0;|r; and dt<t>i\ri etc. are given in terms of ^ via (5.5) 
and analogues thereof. 

Using this result it is straightforward to prove analogous results to those seen in §3 
and §4. We quickly arrive at the following lemma, which is analogous to lemma [2] in the 
case of two spatial dimensions. 

Lemma 5. In m + 1 spatial dimensions, the linearised equations take the form: 



= A/" 2 (*) 



The matrices jVi(fc) and Nzik) take the form 

M(fc)= (^vlk)) ™d M(A) = ( 



o i„ 

Q(fc)-||/c|| 4 I„ -2iU(k) I ' 



where V(k), Q(k), TZ(k) S Mat„(R) are symmetric, tridiagonal and real for k G R' 
The matrix V(k) is given by: 

|/c||+coth(||fc||A/ii)+coth(||fc||A7i ) csch(||fc||A7ii) 

PI ~ " "' " W 

csch(j|fcj|A/ii) ||fc||+coth(j|fc||A/i 2 )+coth(||fcj|A/ii) 



The matrix Q(k) is given by: 

(k-U ) 2 coth(||fc||A/t ) + (fc-C/i) 2 coth(||fc|| Afti) (fc-t/i) 2 csch(||fc||A/ii) 

Pll - " ~ ~ ||fe|| 

(fc-f7i) 2 csch(||fc||Ahi) (k-Ui) 2 coth(||fc||A/u) + (fc-E/ 2 ) 2 coth(||fc|| Ah 2 ) 
Pll \\k\\ 



The matrix lZ(k) is given by: 



(k-U ) coth(\\k\\Ahg) + (k-U 1 ) coth(||fc|| Ahi) (fc- t/i)csch( || fc|| Afai ) 

\\k\\ - \\k\\ 

(fc-t7i)csch(||fc||Afei) (k-Ux) cot!i(|| fc]| Ah 1 ) + (k-U 2 ) coth(||fc|| Ah 2 

Pll 11*11 



It is clear form of the matrices V(k), Q(k) and TZ(k) that the results of §4 & 5 will 
follow in a similar fashion, and as such we simply state the following. 

Theorem 4 (Well-posedness in higher dimensions). Let the n elastic plates be distributed 
horiztonally throughout the channel Q, C R m+1 . If the initial amplitudes i]q and velocities 
dtrja belong to H^ x (R. m ), with s > 2, then the Cauchy problem is locally well-posed in 
L^ x (TL m ). //, in addition, there is zero mean flow, so that that Ui = for i = 0, . . . ,n, 
the problem is corresponding Cauchy problem is globally well-posed. 
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Corollary 1. The linear problem inm+l spatial dimensions is asymptotically stable if 
and only if there is no mean flow, i.e. Ui — for i = 0, . . . , n. 

These results are very much similar to that seen in the two dimensional case. The 
one difference is the conditions needed on the initial data to ensure well-posedness. For 
example, in the physically interesting case of three dimensions, we only need the initial 
data to be continuous and drop the need for differentiablity. This is a consequence of 
the well-known embedding: 

^(R^C^R 2 ). 

In higher dimensions, the conditions are even weaker: it is enough for the initial data to 
belong to L^^ m ^ 4 \ll m ) for m > 4. The case m = 4 is special, and in this instance 
the initial data need be in L^^R 4 ) for each p € [1, oo) to ensure well-posedness. 

The fact that the linear problem is unstable for non-zero fluid speeds is no great 
surprise. Again, the fluid carries vertical momentum between the plates, which will be 
reflected back from the rigid boundary of the channel f2. If the fluid has non-zero velocity, 
and the initial data is non-zero, then the fluid can increase the vertical momentum of 
the plates, and this process is able to continue, allowing instabilities to develop. 



5.4 Periodic Motion 

Problems of physical interest will always involve finite spatial domains, so in certain 
circumstances treatment of the problem on the infinite line is inappropriate. For this 
reason, it is important to be able the analyse the case in which the plates are of finite 
length - one way of achieving this is by considering periodic motion. 

It is well known that for a bounded flows (i.e. flows on compact domains), there is 
in general only a countable infinity of normal modes. As such, one expects the relevant 
global relations, i.e. the integral equations in Proposition [4] etc, to be valid for k in some 
countable subset of R, and this is indeed the case. 

For simplicity assume all functions have period 2ir. In the non-local formulation, after 
constructing a divergence free expression we integrated it over the entire domain f2; (see 



equation (2.3 1 and those thereafter). In the periodic case, the appropriate adjustment 
is to integrate over one cell in f2 i; of width 2tt. In the case of the infinite line, the 
contribution from the vertical components of d£li vanished due to our assumptions on 
the decay of the fields. In the periodic case, we only have the fact that all the functions 
agree on either side of the cell, so to eliminate the corresponding integrals, we must 
choose k so that: 

exp(— ikx) — exp(— ikx — 27rifc), 

i.e. we choose k 6 Z. The modified version of Proposition [4] would now involve integral 
equations valid for k £ Z. Of course, this is sufficient^ to rebuild the functions T]i and 
£f. Indeed, we know from standard Fourier analysis on Lp er that the functions r\i , ^ 
can be reconstructed from their relevant Fourier coefficients, so the global relations, valid 
for k € Z, are simply equations for these coefficients. In the case of the infinite line, 
the global relations can be treated as equations for the Fourier transforms f)i and 
In fact, the global relations can be viewed as pseudodifferential operators acting on the 
Fourier transforms, but we shall not elaborate on this point in this paper. 



3 With moderate a priori conditions on the relevant function spaces. 
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6 Conclusions 



In this paper we have given a detailed study of the linear problem, involving the motion 
of a collection of elastic plates driven by different mean flows. In particular, we study 
the necessary and sufficient conditions needed for asymptotic stability, and also the 
well-posedness of the associated Cauchy problem. In the latter case, we find sufficient 
conditions on the initial amplitudes and speeds of the plates so that the solution of the 
problem changes continuously with the initial data. 

Our approach involves a non-local formalism, in which the underlying system of non- 
linear equations are converted to a coupled system of 1-parameter integral equatoins. 
This approach is largely motivated by a new, unified approach to boundary value prob- 
lems introduced by A.S Fokas in the late 1990's [5J. 

In the two dimensional problem, we show that it is enough for the initial data to 
be continuously differentiable with decay at infinity, for the problem to remain well- 
posed. We then extended the results so that they are applicable to the problem in 
higher dimensions. Using some standard Sobolcv embedding results, we find that in the 
case of three spatial dimensions, the need for differentiability can be dropped, and it is 
sufficient that the initial data be continuous (up to a set of measure zero). 

Finally, we study the asymptotic stability of the problem: for which configurations 
does the Green's function decay to zero after long times, in a frame moving at arbitrary 
velocity? We find that if there is any mean flow, the problem is asymptotically unstable, 
although the unstable modes grow very slowly if the channel is large, and the fluid veloc- 
ities are small. We offer support of this result on physical grounds, based on momentum 
transfer between the fluid and the plates. In the case in which there is no mean flow, 
we show that the Green's function decays like 0(i -1 / 2 ) for large t, which follows from 
analysis of the spectrum of the infinitesimal generator of the semigroup for the Cauchy 
problem in the Fourier space. We also study the pseudospectrum of this operator, and 
discover that for long wavelength modes, the operator becomes highly non-normal. This 
indicates that whilst theoretically the configuration should remain stable for zero mean 
flow, and instabilities should only develop after long time periods in the case of slow 
flow, in practice it is likely that these instabilities will develop at a much faster rate. 
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